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SUMMARY 

According to general relativity, a clock experiencing a shift in the gravitational potential AU 
will measure a frequency change given by A/// « AUjt?. The best clocks are optical clocks. 
After about 7 hours of integration they reach stabilities of A/// ~ 10“^®, and can be used 
to detect changes in the gravitational potential that correspond to vertical displacements of the 
centimetre level. At this level of performance, ground-based atomic clock networks emerge as 
a tool that is complementary to existing technology for monitoring a wide range of geophysical 
processes by directly measuring changes in the gravitational potential. Vertical changes of the 
clock’s position due to magmatic, post-seismic or tidal deformations can result in measurable 
variations in the clock tick rate. We illustrate the geopotential change arising due to an inflat¬ 
ing magma chamber using the Mogi model, and apply it to the Etna volcano. Its effect on an 
observer on the Earth’s surface can be divided into two different terms: one purely due to uplift 
(free-air gradient) and one due to the redistribution of matter. Thus, with the centimetre-level 
precision of current clocks it is already possible to monitor volcanoes. The matter redistribution 
term is estimated to be 3 orders of magnitude smaller than the uplift term. Additionally, clocks 
can be compared over distances of thousands of kilometres over short periods of time, which 
improves our ability to monitor periodic effects with long-wavelength like the solid Earth tide. 

Key words: Atomic clocks; Surface deformation; geopotential measurements, Mogi model; 
Solid Earth tides 


1 INTRODUCTION 

Vertical deformation transients are key to characterising many ge¬ 
ological processes such as magmatic or tectonic deformation (Fig. 
1). Many of these processes have timescales from hours to years 
which are difficult to measure with existing instruments. Atomic 
clocks provide a new tool to resolve vertical displacement, with a 
current precision of about 1 cm in equivalent height after an inte¬ 
gration time of 7 hours (Hinkley et al. 2013; Bloom et al. 2014; 
Nicholson et al. 2015). 

In the past, we argued (Bondarescu et al. 2012) that clocks 
provide the most direct local measurement of the geoid, which 
is the equipotential surface (constant clock tick rate) that extends 
the mean sea level to continents. Since a clock network is ground- 
based, it can provide variable spatial resolution and can be used to 
calibrate and add detail to satellite maps, which suffer from alias¬ 


ing (errors due to effects faster than the sampling rate) and from the 
attenuation of the gravitational field at the location of the satellite. 

In this paper, we consider dynamic sources that cause both 
vertical displacement and underground mass redistribution which 
produce changes in the local geopotential. Geopotential differ¬ 
ences AU are directly measured by the changes in clock tick rate 
A/// ~ AUjc?, where c is the speed of light. To be useful, a 
clock must always be compared to a reference clock, which could 
be nearby or thousands of kilometres away. Clocks are connected 
via ultra-precise fiber links that are capable of disseminating their 
frequency signals over thousands of kilometers with a stability be¬ 
yond that of the clock (Droste et al. 2013). As a concrete example 
we present the case of the inflation (or deflation) of an underground 
magma chamber, computed analytically (Mogi 1958), and apply it 
to the Etna volcano. We explore whether the magma filling could 
be detected using one or two clocks located on the volcanic edifice. 
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Figure 1. Phenomena that could be monitored with optical clock networks. 
The black line shows the lower sensitivity of today’s best optical clocks; see 
equation (1) and Azto^ay in (7). The vertical dashed lines show current and 
planned clock integration times. 


The primary tools currently used to monitor vertical dis¬ 
placement are InSAR and GPS. Interferometric Synthetic Aperture 
Radar (InSAR) measures millimetre displacements in the line of 
sight of radar satellites over wide areas (e.g. Biirgmann et al. 2000; 
Biggs et al. 2011), hut with limited sampling rates (days to weeks). 
GPS is able to measure vertical displacements of 1 cm over short 
timescales (~ an hour) only when the displacement is very local¬ 
ized in the network and/or the frequency of motion is different from 
the frequency of various artefacts that impact GPS accuracy. After 
surveying areas for more than 10 years, the level of accuracy of 
GPS techniques is close to the millimetre level (Blewitt & Lavallee 
2002; Houlie & Stem 2012) or better, enabling us to better char¬ 
acterize the crustal elastic contrast of plate boundaries (Jolivet et 
al. 2009; Houlie & Romanowicz 2011). Since the primary source 
of noise in GPS measurements is due to signal dispersion through 
the atmosphere, both differential GPS and post-processed GPS data 
perform better if networks are dense (e.g. Khan et al. 2010; Houlie 
et al. 2014) because many artefacts cancel across networks over 
which the ionosphere and troposphere can be assumed to be con¬ 
stant. For timescales of seconds, broadband seismometers can be 
used (Houlie & Montagner 2007), but their bandwidth is unsuited 
for resolving long-term displacements (Boore 2003). 

Unlike the GPS or InSAR measurements, local atomic 
clock measurements are insensitive to atmospheric perturbations, 
and could resolve ground displacement over shorter integration 
timescales (hours to months). Further, clocks in conjunction with 
gravimeters are also able to resolve density changes in the Earth 
cmst that do not, or just partially, lead to uplift or subsidence. In 
the case of a spherical magma chamber, the geopotential term re¬ 
sulting from mass redistribution is inversely proportional to the dis¬ 
tance to the source ~ 1/7?, whereas the ground displacement term 
scales with 1/7?^. As opposed to this, both terms have the same 
1/7?^ scaling in gravity surveys. Comparing the measured grav¬ 
ity change to the uplift, Ag/Ah, can reduce this degeneracy to 
model processes a volcano is undergoing before eruption (Rymer 
& Williams-Jones 2000)and also more general processes (see Fig. 
1 and De Linage et al. 2007). 

Even in areas without active seismic or volcanic processes, the 


solid Earth tide has a vertical amplitude that can be as high as 30 cm 
(Agnew 2007) with a semi-diurnal period, whose amplitude can be 
monitored by an atomic clock that is compared with a distant ref¬ 
erence clock. We find that geopotential and gravity measurements 
are sensitive to two different combinations of tidal Love numbers, 
and could be used to calibrate existing measurements of the solid 
Earth tide. 


2 OVERVIEW OE ATOMIC CLOCKS 

According to Einstein’s theory of general relativity, time slows 
down in the vicinity of massive objects. On a neutron star, clocks 
tick at about half their rates on Earth. An observer outside the hori¬ 
zon of a black hole even sees time stopping all together at the hori¬ 
zon. Similarly, clocks that are closer to Earth tick slightly slower 
than clocks that are further away. 

Atomic clocks employ atomic transitions that are classified, 
depending on the transition frequency, as either microwave or op¬ 
tical clocks. Since the clock frequency depends only on a known 
atomic transition and constants of nature, clocks respond identi¬ 
cally to changes in the gravitational potential and do not require 
calibration. This is in contrast to relative gravimeters, which suf¬ 
fer from insh'ument drift and have to be calibrated via comparisons 
with other measurements at the same location. The current defini¬ 
tion of the second is based on a microwave atomic clock. However, 
optical clocks have the potential for higher stability because they 
utilize atomic transitions with resonance line-widths typically 10® 
narrower than microwave transitions. Since the development of the 
femtosecond laser frequency comb, optical clocks have been im¬ 
proving extremely rapidly (Poll et al. 2014). Today’s best clocks 
are optical. They are laboratory devices with frequency uncertainty 

A/// ~ 3 X 10"^®/VT/sec, (1) 

where r is the integration period (Hinkley et al. 2013; Bloom 
et al. 2014; Nicholson et al. 2015), and are likely to continue to 
improve dramatically within the next decade (Poli et al. 2014). A 
transportable optical clock that monitors and compensates for en¬ 
vironmental effects (temperature, pressure, electric and magnetic 
fields) and fits within two cubic meters has recently been built (Poli 
et al. 2014). 

In the future, it is expected that clocks will become sensitive 
to surface displacements at the sub-millimetre level. Possible tech¬ 
niques include nuclear optical transitions (Campbell et al. 2012), 
optical transitions in Erbium (Kozlov et al. 2013), and transitions 
in highly charged ions (Derevianko et al. 2012). A stability of 

^/// ~ CTtomorrow = T / SBC (2) 

should be possible (Hinkley et al. 2013), which would achieve 
A/// ~ 10“^° within one month. 

In order to take full advantage of the improved stability of 
optical clocks, distant clocks have to be compared reliably to the 
10“^® level. This entails a global understanding of vertical dis¬ 
placements, the solid Earth tide, and, overall, of the geoid to the 
1-cm level. Any effects that cause perturbations to this level would 
have to be reliably modelled and understood. Clock comparisons 
via satellites are currently limited by the precision of the communi¬ 
cation link that passes through a potentially turbulent atmosphere. 
The most precise comparisons of distant clocks currently use un¬ 
derground optical fiber links. Optical frequency transfer with stabil¬ 
ity better than the clock has been demonstrated over a two-way dis¬ 
tance of 1840 km (Droste et aZ. 2013), with a fiber-link from Braun¬ 
schweig, Germany to Paris, France. A fiber-link network capable 
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of disseminating ultra-stable frequency signals is being planned 
throughout Europe (e.g. NEAT-FT collahoration, REFIMEVE-I-). 


3 METHODS 


While clocks are sensitive to changes in the gravitational poten¬ 
tial, relative gravimeters see changes in the vertical component of 
the gravitational acceleration, which is a vector (g = — VC/) whose 
amplitude can generally be measured much better than its direction. 
They both provide local measurements of the change in gravity or 
potential relative to a reference point where the gravity and poten¬ 
tial are known accurately. Within a large fiber links network, the 
reference clock could be very far away. 

For a displacement z, the geopotential and gravity changes are 


AU 


GM(b 

K 




2 ,GM^ 


(3) 


where G is the gravitational constant, and and the mass 
and radius of the Earth. The well-known free-air correction Ag for 
gravity is the gradient of the potential change AU. 

A vertical displacement of 2 ; = 1 cm causes changes in the 
geopotential and in gravity of 


AGP 

3 0.l(^) 

V 1 cm/ sec^ 

(4) 

Ag-. 

, 3 X 10-® (^) ^ ~ 3 (^) Mgal. 
VI cm/ sec^ VI cm/ 

(5) 

Thus, the frequency of a clock changes by 



J VI cm/ 

(6) 


The solid line in Fig. 1 is the vertical displacement Aztoday as 
a function of clock integration time r, which is obtained by equat¬ 
ing (1) and (6), while Aztomorrow results from (2) and (6). 


/ j- \-l/2 / j- 1/2 

AZtoday ~ 300 ( ) Cm, Aztomorrow ~ 10 ( J Cm, 

V sec / V sec/ 

(7) 


4 APPLICATIONS 

We first discuss the geopotential and gravity changes caused by an 
inflating magma chamber, using the Mogi model. This is followed 
by a discussion of the measurability of the solid Earth tides. 

4.1 Inflating magma chamber - the Mogi model 

An inflating or deflating magma chamber can be described by the 
so-called “Mogi model” (Mogi 1958): an isolated point pressure 
source in an elastic half-space that undergoes a pressure change. 
The Mogi model is broadly used in the literature (e.g. Ftoulie et 
al. (2006) and Biggs et at (2011)). Recently, its predictions were 
compared to sophisticated simulations, which showed that in many 
situations the discrepancy is quite small (Pascal et al. 2014). 

Clocks and gravimeters lying above the magma chamber will 
be affected by (1) a ground displacement, (2) the change of mass 
density within the chamber, (3) the uplifted rock, and (4) the change 
in the density of the surrounding material. We refer to effects (2), 
(3) and (4) as ‘mass redistribution’ (AU 2 , AU 3 and AU 4 , in Ap¬ 
pendix A). The change of the gravitational potential due to ground 
displacement (1) arises since the observer is shifted to a position 
farther away from the centre of the Earth. A Mogi source centred at 


(0, d) that undergoes a volume change AV deforms the half-space, 
lifting an observer sitting at (r, 0) by 


\w\ 


(l-t 2 ) 


AC d 
TT {r'z -\- d?fG ’ 


(8) 


where v is Poisson’s ratio. Eor crustal rock the typical value is « 
0.25. Using Eq. (3) with z = w, this uplift changes the potential by 


AUi 


GMq _GM(s AC d 

~ "'’AT (r2-fd2)3/2’ 


(9) 


where higher order terms have been neglected. When we add the 
contribution of the mass redistribution terms 


AUm — GpmAV =. (10) 

V -I- 


The total geopotential change is 


AU = AUi -F AUm 


GM® AC d 
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where pm is the magma density. Typically, the change due to mass 
redistribuiton AUm is several orders of magnitude smaller than that 
due to uplift AUi. Measuring the mass redistribution term would 
require that we subtract the uplift term which may be obtained in¬ 
dividually using other techniques like GPS or InSAR. 

Similarly, the change in the gravitational acceleration is given 
by 




Agi + Agm 
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( 12 ) 


Notice that for both the acceleration and the potential the dif¬ 
ferent terms of (2), (3) and (4) mostly cancel each other and only 
a term ~ pm survives (see Appendix A for the complete calcu¬ 
lation). However, this cancellation is a consequence of the elastic 
half-space assumption of the Mogi model and will be less exact 
in more realistic scenarios. On the one hand, for the acceleration 
all terms have the same spatial dependence. On the other hand, for 
the potential they scale differently: while the uplift term scales as 
d/R? like all acceleration terms, the mass redistribution term scales 
as 1/R, where R = is the distance to the source. 

Assuming a clock with stability A/// = (tq/ fr j sec and a 
source with a constant volume change rate AU/Ar, the minimum 
required integration time to measure the uplift at a location (r, d) 


Tl 


aoTtp 

-1 


5 ( 1 - 22 ) 

d ’ 



1 2/3 


(13) 


with g = GM^/Rf, whereas for the mass redistribution 


(JoRp 

(Avy^ - 1 ' 
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(14) 


We focus on the specific example of the Etna volcano. Houlie 
et al. (2006) used GPS data to investigate its underground magma 
system and found that a Mogi source located about 9.5 km below 
the summit with a volume change rate of 30 x 10® m®/yr would 
yield the observed uplift. In Fig. 2, we plot the ground motion and 
the resulting potential change as well as the potential change due 
to mass redistribution. We find that AUm ~ 10“®A(7i, which 
roughly corresponds to A/// ~ AU jt? ~ 10“^®. If we assume 
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Figure 2. Estimate of vertical ground deformation on the Etna volcano over 
the course of one year with the Mogi model. In (a), the altitude profile of 
Etna is shown. The solid curve shows the southern slope of the mountain, 
the dashed line shows the mirror image of the northern slope. For the fol¬ 
lowing plots only the former is considered since the latter gives very similar 
results. Assuming a Mogi source located 9.5 km below the summit with a 
volume change rate of 30 X 10®m^/yr, the vertical and horizontal motion 
over one year is plotted in (b). The change in the gravitational potential due 
to this uplift together with the change due to mass redistribution is shown 
in (c). Notice that the latter is about 3 orders of magnitude smaller. If one 
clock is positioned at the summit, (d) shows the fractional errors on the 
measurements of the source depth (5d/ d and the volume changedAV/I/ as 
a function of the horizontal distance to a second clock. 


constant volume change rate, an optical clock on the summit would 
see the uplift if integrated for about ten days (Eq. (13)), while the 
mass redistribution signal is out of reach. With a clock stability 
(^tomorrow, such an Uplift would be observable within a day, and 
the mass redistribution in five months. 

We also give an example of how to best choose the clock loca¬ 
tions. We assume that the horizontal position of the magma cham¬ 
ber is already known from previous surveys, and that the two mea¬ 
surements we wish to make are the depth of the magma chamber d 
and the volume change /SV (assuming also that v is known). With 
one clock directly above the magma chamber, we find the location 
of a second clock which minimizes the measurement errors on d 
and AV ; of course there must also be a distant reference clock. For 
a flat half-space we find the optimal location to be at a horizontal 
distance r ~ 0.78d (note that the minima is broad; see Appendix 
A5). Assuming Etna to be a half-space with chamber depth 9.5 km, 
the fractional errors there are 5d/d « 18% and SAV/AV « 38%. 
Including the height profile of the mountain. Fig. 2(d) shows 5d/d 


and 5AV/AV as a function of the distance to the summit, assum¬ 
ing a clock sensitive to AU/c? = 10“^®. We find that the opti¬ 
mal location of a second clock is r « 0.78 d, independent of the 
clock’s performance. For the performance considered here, this cor¬ 
responds to &d/d « 14% and SAV/AV ~ 30%. 


4.2 Solid Earth Tides 

Solid Earth tides are the deformation of the Earth by the gravi¬ 
tational fields of external bodies, chiefly the Moon and the Sun. 
In general relativity the Earth is freely falling, so at the centre of 
the Earth the external gravitational force is canceled by the Earth’s 
acceleration toward the external body. Because the external field 
is not uniform, other points experience a position-dependent tidal 
force (Agnew 2007). 

The tide has three effects of concern on the gravitational po¬ 
tential. First is the external potential itself, which can be calculated 
directly from the known mass and position of the external bodies. 
Second is the change in the Earth’s gravitational potential produced 
by the deformation of the Earth. Third is the vertical displacement 
of the surface, which produces a free-air correction to our measure¬ 
ments. The last two effects are proportional to the first, parameter¬ 
ized by the Love numbers kn and hn, respectively (Agnew 2007). 

As shown in Appendix B, clocks (sensitive to geopotential 
changes) and gravimeters (sensitive to the downward component 
of g) each measure a different linear combination of the Love num¬ 
bers. This is because the three effects scale differently with the dis¬ 
tance from the centre of the Earth. It is therefore desirable to com¬ 
bine gravimeter and clock measurements to infer the Love numbers 
with great accuracy. Recall that clock measurements must be made 
between a pair of clocks. Since tidal effects are global, both clocks 
are sensitive to the tidal deformation; one clock cannot simply be 
treated as a reference. To measure the tidal amplitude, it is nec¬ 
essary to compare instmments over distances of the order of half 
the tidal wavelength on timescales shorter than the period of the 
tide. The tidal wavelength is half the circumference of the Earth for 
the dominant tidal mode. Both differential GPS and InSAR have 
short baseline, so tidal effects are small and are mostly subtracted 
by modelling. A network of clocks where each clock pair is sep¬ 
arated by hundreds of kilometers could more holistically monitor 
solid Earth tides. Current clocks provide measurements of the ver¬ 
tical uplift to within a percent of the maximal tidal amplitude on an 
hourly basis. Such measurements could be used to monitor stress 
changes within the crust, and to investigate whether these correlate 
with triggered seismicity. 

Such a network of clocks on the continent scale could accu¬ 
rately measure the tidal Love numbers. The external tidal potential 
can be decomposed into a sum of Legendre polynomials 

Utid{R,a) = - - — I — 1 Pn(cosQ), (15) 

where R is the distance between the external body of mass Mext 
and the Earth’s centre of mass, and a is the angle from the ob¬ 
servation point to the line between the centre of the Earth and the 
external body. Since 77® /77 « 1 /60 for the Moon, and 7^0 /77 « 
1/23000 for the Sun, it is typically sufficient to just consider the 
first few terms of the expansion. By linearity, we treat each n sepa¬ 
rately. 

The potential change measured by a clock, including all three 
effects above, is 

A(/„ = (1 + -/i„)(/"a(77, a) . 


(16) 
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The change in the vertical gravitational acceleration measured by a 
gravimeter is 


Ag„ = -n 



n + 1 
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kn “h — hr 
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Um{R,c^) 

R® 


(17) 


Combining the two measurements with the known <^) we 

can determine the Love numbers 
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(18) 


(19) 


Both Love numbers are believed to be modelled to within a fraction 
of a percent (Yuan & Chao 2012). These models will be throughly 
tested once a global optical clock network becomes available. 


5 CONCLUSIONS 

We have demonstrated the promise of very precise atomic clocks 
for geophysical measurements with two illustrative examples. In 
the inflation or deflation of a spherical magma chamber (the Mogi 
model) clocks are primarily sensitive to the local vertical displace¬ 
ment resulting at the Earth’s surface. Such monitoring of local de¬ 
formations can be done using a reference clock anywhere outside 
the zone where the displacement is significant, typically tens of 
kilometres. However, when monitoring solid Earth tides, the best 
accuracy will be obtained with a clock network spanning the globe. 

Beyond the examples given, it should be possible to use clocks 
in conjunction with gravimeters to monitor dyke intrusion, post- 
seismic deformation, aquifer variations, and other effects causing 
vertical displacement or subsurface mass change. In contrast to 
GPS and InSAR measurements, ground clock measurements are 
insensitive to the turbulence in the atmosphere. 

In the future, optical clocks are expected to become part of 
the global ground-clock network that is used for telecommunica¬ 
tions increasing the precision with which we can monitor time. 
Ground clocks can thus be combined with existing instrumenta¬ 
tions (GPS, InSAR, gravimeters) to track underground mass redis¬ 
tribution through its effects on the geopotential. Portable optical 
clocks have already been developed (Poll et al. 2014), and could 
monitor changes in the geopotential of the Earth across fault lines, 
in areas with active volcanoes, and for surveying. 
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Figure Al. Deformation due to inflating cavity. 


APPENDIX A: T HE MOGI MODEL IN DETAIL 


A Mogi source is an isolated pressure point source embedded in an 
elastic half-space. A pressure change AP deforms the surround¬ 
ing half-space and uplifts an observer standing on the surface. This 
ground displacement, discussed in Sec. 4.1, affects the ticking rate 
of a clock since its position in the Earth’s gravitational field is 
slightly shifted. Besides that, there are three additional effects af¬ 
fecting both the gravitational potential and acceleration, which are 
due to the redistribution of matter. These are discussed in Sec. A2, 
A3 and A4. 

We choose coordinates such that the elastic half-space ex¬ 
tends from minus to plus infinity in both the x and y directions, 
and from 0 to -|-oo in the 2 direction. The source is located at 
(0, 0, d) and an observer positioned on the surface has coordinates 
{x, y, 0). By symmetry, we typically use coordinates (r, z), where 
r = \/x'^ -\- y'^ is the distance from the observer to the place on 
the surface directly above the cavity. A scheme is given in Fig. Al. 

Even though a point source is considered, the pressure change 
can be interpreted as the volume change AV of a finite size cavity. 
This approximation is valid as long as the radius of the cavity is 
much smaller than the depth d. Such a volume change can, for ex¬ 
ample, be caused by the inflation or deflation of a magma chamber. 

The radius of the cavity changes by (Hagiwara 1977) 


Aa 


aAP 


(A.1) 


where y, is the shear modulus, having the same units as pressure 
(note that some authors use G for the shear modulus). Typical val¬ 
ues for rock are ^ = 10 x 10® Pa to 30 x 10® Pa. Thus, the volume 
changes by 

AF = ^{a -f Aa)® - ya® = 47ra®Aa -f C>(Aa®). (A.2) 

Neglecting terms of order 0{Aa^), it can be written as 


AV" 


Tra® AP 


3V"AP 


(A.3) 


The change in pressure within the cavity induces a volume change 
that affects the surrounding elastic medium. An observer on the 
surface at (r, 2 = 0) is displaced by 



= a®AP 




(l-u) 


AV" 

TT 



(A.4) 
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Figure A2. Integration times for current clocks. Considering a clock with 
ctq = 3 X 10“^®, the required integration time to resolve uplift (Eq. (13)) 
is shown in (a) as a function of distance R and (constant) volume change 
rate AV/ At. Here, it is assumed that the clock is placed above the magma 
chamber (r = 0). Analogous, the integration time for seeing the redistribu¬ 
tion of mass (Eq. (14)) is shown in (b). In contrast to the uplift where the 
result depends on both r and d, the integration time for the mass redistribu¬ 
tion term depends on the total distance R only. 


with u being the Poisson ratio and R = yjr'^ + being the dis¬ 
tance between the cavity and the observer. Thus, a point on the sur¬ 
face that was originally at (r, 0) gets shifted to (r -h u, tu). There¬ 
fore, an observer on the Earth surface is vertically uplifted by 

l^(")l = (l-")^ (r^ + rf^)3/2- (A.5) 

In addition, the density of the elastic body changes at each 
point by (Hagiwara 1977) 

pa® AP r® - 2{z + df 


Ap{r,z) = 


= {l- 2 v)p- 


where A is Lame’s constant 


X + y [r'^ + [z + 

AV^ r®-2(z-Vd)® 


TT [r^ -\- {z + d)2]5/2 ’ 


_ 2 yv <'=1 

X = - = y. 

l-2u ^ 


(A.6) 


(A.7) 


Below, we calculate the changes in the clock and gravime¬ 
ter measurements due to ground displacement and the three effects 
coming from mass redistribution. Summing the results of Appendix 
A2, A3 and A4 we find 

AP„ = AU 2 + AU 3 + AUi = -Gp™AV" 


\/r^~+~d? 


(A.8) 
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with the expansion 


g{r,z) 


GMq 2GM0 GMq . 2 


3z^) . 


(A. 12) 


The linear term in the expansion is known as the free-air cor¬ 
rection; we can neglect the higher orders. We define the gravita¬ 
tional acceleration at the surface as, g = GM^ / Rf « 9.81m/sec^. 
For a Mogi source z = w (see Eq. A.4). The potential and gravity 
changes are 


AUi 

Agi 


= -g z = g{l - v) -. 

TT 

= 25 ^ = - 2 p(l - u) 

A® 


d 

(r2 -f d2)'3/2'’ 

Al/ d 

TT 7?® (r2 -(- ((2)3/2 


(A. 13) 
(A. 14) 


A2 Direct signal 

If there is inflow or outflow of magma, crustal rock of the vol¬ 
ume AV will be replaced by magma, or vice versa. This change 
in density gives a direct signal, due to the mass change AM — 
AV{pm — p), where pm. and p are the densities of magma and 
crustal rock, respectively. The resulting change in the gravitational 
potential from the presence of the density anomaly only is 

AU,.-gfi.-G(p,,-p)AV-^^. (A.15) 

and the change in the gravitational acceleration is 

^ ^ - P^^^ id^+%/2 - (A.16) 


Figure A3. Integration times for future clocks. Analogous to Fig. A2, the 
integration times to resolve uplift and mass redistribution are shown for 
future clocks with (tq = 1 X 10“®®". 


The required integration timescales to resolve uplift (n) and mass 
redistribution (Tm) as a function of the distance to the magma 
chamber and the volume change rate AV /At is given by Eq. (13) 
and (14), respectively. Eor a current clock with ao = 3 x 10“®®, 
the respective integration times are shown in Eig. A2(a) and (b). 
Resolving A[/m requires either long integration timescales of the 
order of a year or more or better clocks. However, for different ge¬ 
ometries (e.g. flatter magma chambers) or more realistic models 
the mass redistribution term could be more significant. Eor future 
clocks with ao — 1 X 10“®®^, the analogous plot is given in Eig. 
A3. 


Al Ground displacement 


Treating the Earth as a perfect sphere, its gravitational potential is 


U(r,z) 


GM(s 


(A.9) 


where M® is the mass and 77® is the radius of the Earth. 

We expand the potential in a Taylor series around (r, z) = 
(0, 0) and neglect all terms of order 0 {r^), 0{z^) and higher to 
obtain 


U{r,z) 


GAT® 

R@ 


GM^ GM^ , 2 


20^) 


(A. 10) 


The gravitational acceleration is given by 


gf, z) 


GM^ 

r2 -F (77® - z)2 ’ 


(A.11) 


A3 The potential of the uplifted rock 


Before the inflation of the cavity, the surface of the half-space is 
flat. But after the inflation there will be a hat of material peaking 
above the location of the cavity; this is where the clock and the 
gravimeter are located. Obviously, this hat of material affects the 
outcome of a measurement. As a simplification, most authors re¬ 
place the hat by an infinite disc of height do. While this gives the 
correct first order result for the gravitational acceleration, the ex¬ 
pression for the gravitational potential diverges. 

In this section we discuss the gravitational potential and accel¬ 
eration of this additional hat without making the disc approxima¬ 
tion. This allows us to calculate the effect on clock measurements 
and to compute higher order corrections to the gravitational accel¬ 
eration. 

We denote the uplift directly above the cavity by 

do = |w(r = 0)1 = (1 - iz) —^ . (A.17) 

TT 

We shift the coordinates slightly to put the observer at z = 0. Thus, 
z = do now corresponds to the surface of the ground before the 
uplift. The z coordinate of the new surface at a given distance r is 


z(r) = do — |w(r)| 


= do 




-3/2‘ 


This can be inverted to express r as a function of 2 : 


( 5 ) 


, ( , N -W 


-A) = I 1- X 

do 


- 1 . 


(A. 18) 


(A. 19) 


The z index is used to emphasize that is the distance at which 
the surface has altitude z. 
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In general, the gravitational acceleration and the gravitational 
potential of an observer located at (r, z) = (0,0) due to the density 
distribution p(r, z) are 

n(r z)z 

53 = 2nG I dz, (A.20) 

U3 = -2TvGf f -J^dl^rdrdz. (A.21) 

J —oo Jo v 7* Z 

Assuming constant density, the acceleration Ago and the potential 
AUs due the uplifted mass are 

fdo rr^ 

Ag3=2nGpJ^ (A.22) 

fdo r^z „ 

AUo = —2 'kGp / / zdrdz. (A.23) 

Jo Jo Vr^ + 2^ 

Performing the integrals over r yields 

A» . - l) (A.24, 

AU 3 = —2'itGp j ^\/r^ + z^ — zj dz, (A.25) 


and defining a new variable (^ = a/do, we obtain 



= - 2 nGpdodj^J + dC- 

(A.27) 


Because do/d ^ 1, we can expand around do/d = 0 and drop all 
terms that are of order 0 {d^/d?), then we can integrate and find 




2 'KGpdo 




do 

d 


AUo 


— 2 'KGpdod 1 


1 do 
2 ~d 


(A.28) 


To combine this with the other mass redistribution terms, we need 
the result as a function of r. In the infinite disk model in the litera¬ 
ture, the gravitational acceleration at r is taken to be that produced 
by a disk with a thickness equal to |ui(r)|, the local value of the 
uplift. 


Ap3 


27rGp|w(r)| 


(A.29) 


Integrating with respect to d we obtain the corresponding potential 
AUo = -2Gp{l - y)AV^^:=. (A.30) 

V -I- d^ 


A4 Gravitational potential due to density change 


In response to the inflating magma chamber, the density through¬ 
out the surrounding rock changes by Ap{r, z), given by (A.6). The 
resulting change in the gravitational acceleration 

Ap4(r) = -G(l-2^)pAV^-^^^. (A.31) 

was calculated by (Hagiwara 1977). Following a similar approach, 
we determine the change in the gravitational potential due to this 
density variation. This may be useful for the reader since Hagi- 
wara’s paper was written in Japanese. 

First, a new function K is defined 


Kix,y,z) = {x^+y^+ {z + df) , 
which allows to rewrite the density change (A.6) as 


Ap(a;, y, z) = {1 - 2v)p— 

TT 


K{x,y,z) + {z + d) 


(A.32) 


dK{x,y,z) 


dd 

(A.33) 


The gravitational potential due to this change in density is 

GAp{x', y', z')dx' dy' dz' 


AU 4 ,{x,y,z) = - 


r3 ^{x — x'Y -I- (p — y'Y + {z — z’Y 
d 


1 (z -f d) 


dd 


^K{x,y,z,z')dz , 

(A.34) 


where we defined 


^K{x,y,z,z) = 

Jr 


K{x', y', z')dx' dy' 


R2 \/{x — X'Y + {zy — yiy + {z — Z'Y ’ 

(A.35) 


K = -G{l-2v)p—. 

TT 


(A.36) 


Further, we assume that before the inflation of the magma chamber 
the density p is constant. 

The 2D Fourier transform, assuming 2 ; = 0, is 




c(k^,ky) = ^ [ ^K{x,y,0,z')e dp 

277 Jr2 

f K{x', y', dy' 

Jr3 


-f 

277 Jll2 


^ — i(kj:U+kyv) 


du dv 


= 271- 


= 27r 


[27r Jg2 _(_ ^2 _|_ ^'2 

(,-W + d\y/W^ g-|2'lx/fej+fej 

TTd ^ki + ki 

g-{2z'-|-d)x/fej+fcj 


( 2 ' -I- d)^kl + kl 


(A.37) 


We substituted u = x—x' and v = y—y' and used dx dy dx' dy' — 
du dv dx' dy'. Further, we used the Fourier integrals 


g-KkzzX + kyv) g-NIx/'^i+feJ 

o / —T3— 2 ^dxdy= — , (A.38) 

277 J^2 (a;2 +y^ + z^yG -f kl 

-i(ky:X+kyy) -|z|.^fc2+fe2 

dxdy = -, (A.39) 


-f 

277 J ^2 


277 Jk 2 (a;3 + y ^ + 22)3/2 - ^ 

and we used that the density vanishes for z' < 0. Taking the Fourier 



Atomic clocks as tools to monitor vertical surface motion 9 


transform of A [/4 


Af /4 Ikx^ — 


27r . 


K2 


AU4,ix,y)e 


— i{kxX + kyy) 


dx dy 


1 + (2: d) 


- [ ^K{x,y)e-^^’^‘^^+’^y^'>dxdy] dz' 
StT J^2 J 

^*K{k:c,ky)dz , 


l + {z' + d)- 


and using the expression for above, we obtain 


(A.40) 


/ 


AU4{kx,ky) = 2'kk 
= —2'KK 


l + (2' + d)^ 


(2z'+d) 


dd\ (^z'+ d)fkl + kl 

g-(2z' + d)i/fci+feg^^/ 


dz 


= —KK - , 

■fk^+ky 

{AM) 

where we used that outside the half-space (i.e. for z' < 0) the 
density (contained in k) vanishes. The inverse Fourier transform is 

AUA{x,y) = ^ [ AU4{k^,ky)e^^'°'^'^'^^^^dkxdky 
2^ Jr2 


-d-^kl+k^ 

2 7r2 ^kl + 


^ Jm 


^^k^^+'^yy)dk^dky, 


(A.42) 


and introducing polar coordinates k^ = k cos 9,ky = k sin 9 (such 
that dkxdky = kdkdff) and orientating the coordinate axes such 
that {x, y) = (r, 0), the integral can be written as 


AU 4 {x,y) = --K 
1 

=- K 

2 


foo p27v —dk 


0 ^0 
poo 

I —dk 


e“’'“=%0Mfc (A-43) 


Flere, we used the definition of the Bessel functions and the integral 
6.751(3) in (Gradshteyn & Ryzhik 2007). Finally, we find 

1 


AU 4 {r) = G(1 - 2v)pAV 


fr'^ + d? 


(A.44) 


A5 Optimizing the placement of a second clock 


We now consider the question of the optimal location for measur¬ 
ing the depth and volume change of a magma chamber. We will 
assume that one clock is directly over the magma chamber. What is 
the optimal location of a second clock? Recall that there must also 
be a distant reference clock. For this example, we assume that the 
horizontal location of the magma chamber is known from a previ¬ 
ous survey. If this were not known, it could also be measured, but 
more clock locations would be required. 

The shift in the gravitational potential at clock i can be written 


AUfr) = CqAV 


d 

(rf -b d2)3/2 ’ 


_ GM® 1 — V 

^0 = - 

i?® K 


(A.45) 

(A.46) 


where d is the depth of the magma chamber, Vi the horizontal dis¬ 
tance of clock i from the magma chamber. We assume that iz is 
known, thus Go is a constant. 

We measure two quantities, the differences between each of 
the clocks and the reference clock. These are 

AUi = CoAV^, (A.47) 

d^ 

AU, = CoAV ^^^ (A.48) 

The measurements of each of these quantities is affected by 
clock noise. The three clocks will have uncorrelated noise, but the 
AUi will be correlated because they each depend on the reference 
clock. We write the variance of the two clocks near the magma 
chamber as af, and the variance of the reference clock as a\. 
Then the variances and covariance of A[/i and AU 2 are 


Var(A(7i) ^ erf + af 

(A.49) 

Var(A(72) = erf + erf 

(A.50) 

Cov(A(7i, A(72) = . 

(A.51) 


The noise variance in the measurement of d is 

(A.52) 

and the same holds replacing d with AV. It is easiest to calculate 
the partial derivatives by calculating first the partial derivatives of 
the AUi with respect to d and a, and inverting. They are 

dAV _ d^ir'^ - 2d?) 
dAUi ^ 3Gor2 ’ 
dAV _ 2{r‘^ + d?)^'^ 
dAth ~ 3Godr2 ’ 
dd _ d^{r^ + d^) 

dAUi ' “ 3GoAl/r2 ’ 
dd _ (r^ + d^ f!'^ 
dAth ^ 3GoA’l/r2 ’ 

We find the optimal location for the second clock by minimiz¬ 
ing the variance of either AV or d, as a function of r, the horizontal 
distance to the second clock. Fortunately, we find that d and AV 
are optimized for similar values of r, and the minima are fairly 
broad. A reasonable compromise, in the case where the reference 


clock has the same performance as the other clocks, is r ? 
With this clock placement, we have 

s 0.78d. 

/ d^ " 
Var(Al/) « 4.6 — 

\ <^02 

2 

\ 2 

1 

(A.57) 

/ d^ " 

Var(d) « 2.2 

2 

\ 2 
ctm- 

(A.58) 


We can summarize this more succinctly in terms of fractional 
errors. We define 5h = aulg, which is the standard devia¬ 
tion of the measurement of a single clock, in terms of equiva¬ 
lent height. The maximum value of the uplift (at the summit) is 
h = CoAV/{gd?). The standard deviation of one clock is then 
Cm = {C(,AV/d?)5h/h. Given this, the fractional accuracies of 


(A.53) 

(A.54) 

(A.55) 

(A.56) 
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our measurements are 


5/W 

AV" 

5d 

It 



(A.59) 

2 A 

h 

(A.60) 


A frequency stability of 10“^® corresponds to 1 cm, which is 
achieved with current clocks after an integration of 7 hours. If each 
clock has this performance and the maximum uplift is 10 cm, this 
is a fractional accuracy of 10%, giving a fractional accuracy on d 
of 22% and on AV^ of 46%. If instead, we integrate for a month, 
the fractional accuracies improve by an order of magnitude leading 
to fractional accuracies of 2.2% and 4.6%, respectively. 

In the case of a reference clock that is much better than the 
others (gr gm), the numbers become r « 0.87d, 5AV/AV « 
3.7 5h/h, and 6d/d ~ 2.06h/h. 


APPENDIX B: SOLID EARTH TIDES 


Throughout this section we follow (Agnew 2007) and we use stan¬ 
dard spherical coordinates. The gravitational potential of an exter¬ 
nal body can be written as a multipole expansion 


Ucxt 


P 

R 


n=0 


(B.l) 


Here, p is the distance between the external body and the location 
of the clock, R is the distance between the body and the Earth’s 
centre of mass, r is the distance between the Earth’s centre and the 
clock and a is the angle between the two vectors pointing from 
the centre of the Earth to the external mass and the clock. The 
denote the Legendre polynomials. In the sum, the n — 0 term can 
be neglected since it is constant and therefore does not contribute to 
the force since it will drop out when taking the gradient. The n = 1 
term (17ext,n=i = — GMexi/ R^r cos a) causes the orbital accelera¬ 
tion and therefore, by the definition of tides, the tidal potential is 


f/tid(r, q) = P„(cosa). (B.2) 

n=2 


By Ptld we will denote the n-th order component of the expansion. 

The tidal deformation is mainly caused hy the Moon and the 
Sun; the other planets cause very minor effects. Since typically r/R 
is small (for the Moon: r/R ~ 1/60, Sun: r/R ~ 1/23000), it is 
sufficient to just consider the first few orders of the expansion. 

If the Earth was a completely rigid spherical body with no 
external mass acting tidally, the gravitational acceleration at the 
surface would be a constant g = GM^/R/g. Changing the ra¬ 
dial distance by a small amount Ar would change g by Ag = 
— 2GMQ/R^Ar. Thus, going farther away (Ar > 0) weakens 
the gravitational acceleration. 

First, we consider just the effect of an external body, letting the 
Earth be perfectly rigid. This corresponds to the Love numbers (in¬ 
troduced below) being zero. The value of the tidal potential on the 
surface is ?7tid(P®, cx). This induces a change in the gravitational 
acceleration at the surface of 

Ag{a) = - 

dr 

_ _ q) 


Just taking the first (n = 2) term, this is 


Ag{a) 


R® 

P®(3cos^a - 1). 


(B.4) 


Thus, there is an outward pull, Ag > 0, at the sides of the Earth 
facing (a = 0) and opposing {a = n) the external body. Halfway 
between, there is an inward pull (a — 7r/2, 37r/2). Since the Earth 
is not a rigid body, it is deformed by tidal forces. Its tidal response 
is quantified by the Love numbers and hn ; there is also a Love 
number which we do not use here. The larger the Love numbers, 
the stronger the deformation. 

The gravitational acceleration arising from the tidal interac¬ 
tion has three components. First, there is the direct effect from the 
external body (B.3). Second, the deformation of the Earth gives rise 

(n) 

of an additional gravitational potential k„U/^ , resulting in an ad¬ 
ditional acceleration. Third, due to the deformation, the altitude of 
the observer changes by —hnU^//'' /g compared to a non-deformed 
Earth. Thus, the effective tidal potential measured by a clock on the 
surface, for a given n, is 


= (1 + fcn - • (B.5) 


The potential of the deformed Earth is fcn(7(n) • It can be writ¬ 
ten as an expansion around its value (7tid(P®, a) at the undistorted 
surface. Since the gravitational field satisfies Poisson’s equation, it 
is a harmonic function outside the Earth and therefore, we can write 
it as an expansion of the form 


[/Earthdef(r, «) = I] ( ^ ) Fr,{a), (B.6) 

n=0 ' ^ 


where the functions F„{a) have to be chosen such that they match 
Utid{R^, a) at the Earth’s surface, giving 


[/Earth def(t*, Ct) — 


GMe, 

R 


E(%) (^) 

(B.7) 


n + 1 


The tidal acceleration is now given by 

d{Utid + fcn [/Earth def) 


dr 


GM„ 

R 


n=2 

X P„(cosa) 
which can be written as 


f;^„._Ma + t,)(^)' 


-E" 1- 


n -I- 1. 


Udd{R®,(^) 

R(n 


(B.8) 


(B.9) 


The change in the gravitational acceleration due to a 
displacement Ar close to the surface of the Earth is 

2GM® ^ _ ^GM^ Ar _ ^ Ar 

Rl Rl 77 ® - ^^ 77 ®- 

The final effective tidal acceleration on the surface is 

A5t^®L = (n - (n + l)fc„ + 2h„) 


small radial 


(B.IO) 


(B.ll) 


From (B.5) and (B.ll), we see a key feature of clocks relative 
to gravimeters. Gravimeters are relatively more sensitive to effects 
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at higher n, i.e. shorter wavelengths, than are clocks. This is be¬ 
cause they measure a spatial derivative of the potential rather than 
the potential itself. In contrast, the response to uplift only depends 
on motion in the entire gravitational field of the Earth, and so the 
coefficients of hn do not depend on n. 

Note that clocks are best at constraining the lower order mul¬ 
tipoles. Instruments that measure higher order derivatives of the 
geopotential (e.g. gradiometers) are more sensitive to the higher or¬ 
der multipoles. Clocks will only be useful for measuring the n — 2 
and possibly n = 3 tides. However, a detailed analysis of the mul¬ 
tipole structure based on clock data is beyond the subject of this 
paper. 



